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Abstract 

We study both the spin-average and spin-dependent structure functions of 
the Uthium isotopes, ^~^^Li, which could be measured at RIKEN and other 
nuclear facilities in the future. It is found that the light-cone momentum 
distribution of the valence neutron in the halo of ^^Li is very sharp and sym- 
metric around y = 1, because of the weak binding. This implies that such 
neutrons are in an environment very close to that of a free neutron. The EMC 
ratios for Li isotopes are then calculated. Furthermore we investigate a new 
ratio, Rd, of the difference between the Li structure functions of mass number 
A (-^2^) and A — 1 {F,^~^) to the difference between the structure functions 
of the deuteron and free proton. We study the possibility of extracting the 
neutron structure function from data for the nuclear structure functions of 
the Li isotopes. In particular, the ratios i?„ for ^ = 9 and 11 present an 
attractive possibility for extracting the free neutron structure function. Next 
we calculate the spin-dependent structure functions of ^'^'^^Li isotopes, which 
have spin of 3/2. For such nuclei the spin structure function is given in terms 

3 /2 

of the multipole spin structure functions; ( gi, which is analogous to the 
usual spin structure function, gi, for a target with spin 1/2, and a new one, 
^3^51) which first arises for a target with spin 3/2. The effect of the nuclear 
binding and Fermi motion on ( 171 is about 10% in the region x < 0.7, but 

3 /2 

it becomes quite important at large x. The spin structure function of 3 gi 
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is negative at small x but it becomes positive in the region 0.2 < x < 0.5. 
However, the magnitude is very small. At large x it is again negative and 
its absolute value becomes large because of the Fermi motion. Finally, we 
discuss the modification of the Gottfried and Bjorken integrals in a nuclear 
medium and point out several candidates for a pair of mirror nuclei to study 
the flavor-nonsinglet quark distributions in nuclei. 

PACS: 25.30.Mr; 13.60.Hb; 24.85.+p 

Keywords: Nuclear structure functions; Neutron structure functions; Unsta- 
ble nuclei 
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Typeset using REVT^ 



I. INTRODUCTION 



The recent development of secondary radioactive beams with sufficient intensities, which 
are produced through high energy fragmentation of nuclei, has opened a new frontier in 
nuclear physics, i.e., the study of nuclei with extreme ratios of the numbers of neutrons, N, 
to the number of protons, Z. Such nuclei, far from the stability line, have been investigated 
intensively at many nuclear facilities in the world as one of the main subjects of nuclear 
science today . There are various of interesting problems such as weakly bound nuclei, low- 
density nuclear matter, novel nuclear structures of halo and skin, very asymmetric nuclear 
matter, appearance/disappearance of magic numbers in the vicinity of the neutron and 
proton drip line 

In particular, it is noticeable that the halo structure in some light nuclei was discovered 
on and close to the neutron drip line 0. Those nuclei have spatial extensions much larger 
than naive expectations based on their mass numbers. A typical example is ^^Li, which has 
a root mean square (rms) radius of about 3.1 fm. In Table | we summarize the observed 
data of ^"^^Li isotopes The experimental data suggests that nuclei like ^^Li could 

be understood as a core surrounded by one or two more neutrons, which extend to several 
times the nuclear radius - the so-called halo nuclei. This picture is also supported by recent 
experiments with break-up of loosely bound neutron-rich projectiles on light targets [0. 

The /3-unstable nuclei in the vicinity of the neutron drip line have two distinctive features. 
The ffist feature is a large asymmetry in Z and N. For stable nuclei with Z > 3, due to the 
Coulomb interaction and the nuclear saturation property, (i) the ratio N/Z is 1.0 ~ 1.6, (ii) 
the observed central density is ~ 0.17 fm~^ except for very light nuclei, and (iii) the density 
distributions of protons and neutrons are similar and the thickness of the nuclear surface is 
almost unchanged for all nuclei. For /3-unstable nuclei close to the neutron drip line, the 
N/Z ratio is about 2 ~ 4 and the distributions of protons and neutrons are considerably 
decoupled (i.e., neutron halo and neutron skin). 

The second feature is a small separation energy of the valence neutron. For light stable 
nuclei, where N Z, the separation energy of the last nucleon is usually several MeV. On 
the other hand, with an increasing number of neutrons the separation energy of the last 
neutron, S'„, varies from several MeV to less than 1 MeV at the neutron drip line, while 
the separation energy of a proton tends to increase and sometimes exceeds 10 MeV near the 
neutron drip line. This difference between the separation energy of a neutron and that of 
a proton implies a difference in the Fermi energies of neutrons and protons, and it provides 
an understanding of the mechanism of formation of neutron halos and skins as well as the 
instability of /3-unstable nuclei. 

The partonic distribution functions in such nuclei are of considerable interest in view of 
these features 0. The distributions of quarks in a nucleus differ significantly from those 
in the free nucleon - the so-called European Muon Collaboration (EMC) effect [§,0. Since 
the early 80's a large amount of data on flavor-singlet nuclear structure functions, i.e., 
N c:^ Z, has been collected. However, relatively less attention has been paid to the flavor- 
nonsinglet structure functions of a nucleus. Measurements of nonsinglet nuclear structure 
functions in deep inelastic scattering (DIS) could shed light on many phenomena involving 



nonperturbative QCD in nuclei, such as SU(3) symmetry breaking, flavor asymmetry |T0 
and so on. 
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Recently we have studied the nonsinglet structure function of the hghtest mirror nuclei 
(^He and ^H) 0,0 and pointed out that the nonsinglet structure function in those nuclei 
is enhanced at small x by nuclear shadowing. Furthermore, we suggested that the Gottfried 
integral is generally divergent because of charge symmetry breaking. We have also investi- 
gated the effect of the medium modification of pion fields on the flavor-nonsinglet structure 



function of the lightest mirror nuclei |]12[ . It is intriguing to study such nonsinglet structure 
functions using a pair of unstable nuclei |^ because such nuclei could give rise to a large 
asymmetry in isospin space - recall the first feature for unstable nuclei discussed above. The 
current rapid development of radioactive nuclear beam techniques at many nuclear facilities 
in the world will provide good opportunities for advancement to the study of nonsinglet 
structure functions in the near future . 



Utilizing the second feature of neutron-rich nuclei far from the stability line, i.e., a 
wide range of the neutron separation energy, we could investigate the binding (off-mass 
shell \Ti\ and shell structure dependences) and Fermi motion effects on nuclear structure 



functions 0. Furthermore, the neutron structure function, F2, could be extracted from 
nuclear structure functions of neutron-rich nuclei, by virtue of the weakness of the neutron 
separation energy |^. 

It is very important to determine the neutron structure function from nuclear experimen- 
tal data. The proton structure function, , which at intermediate and large-x measures a 
charge-squared weighted combination of the valence u and d distributions, depends mostly 
on the u quark distributions owing to the larger weight of the u quark charge. It is therefore 
necessary to know the neutron structure function as precisely as possible to determine the 
individual isospin distributions separately. 



The problem of extracting neutron structure functions from nuclear data is rather old [jT5 
Because of the absence of a free neutron target the deuteron is traditionally used as an 
effective neutron target: ^ F2 — F2 {F2 being the deuteron structure function). While 
this approximation may be valid for 0.1 < x < 0.4, it breaks down dramatically for large x 



because of the Fermi motion. The Adelaide group has recently reanalyzed the deuteron 
structure function to determine F^, and they showed that the nuclear binding and Fermi 
motion effects were significantly larger than those in earlier analyses. In particular, ignoring 
the effect of nuclear binding yields an error of up to 50% in the ratio F^/Ff extracted at 
X ~ 0.75. 

Recently, the A = 3 system (^He and '^H) is also being seriously considered as an al- 
ternative way to extract the neutron structure function |T^,|TB[. Using realistic Faddeev 
wave functions, it was shown that F^ could be extracted from the ratio of those structure 
functions, where the nuclear effects are cancel each other so that the ratio is within 2% of 
unity for x < 0.85. In the present paper we discuss another way to extract the neutron 
structure function using the radioactive Li isotopes. This is based on the novel idea that the 
valence neutrons in the /5-unstable nuclei are so loosely bound that they may be regarded 
as essentially as free neutrons. 

For the last decade polarized DIS experiments have yielded a number of important re- 
sults concerning the spin structure of the nucleon. The measurements of the proton spin 
structure function, g^, by EMC [0, combined with flavor-nonsinglet matrix elements from 
weak decays, gave signiflcant information on the singlet axial charge of the proton. The 
small size of this resulted in the so-called proton spin crisis, and it has been reanalyzed 
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by many groups [0. As an extra source of information, it is important to measure the 



neutron spin structure function, g^, for testing the fundamental Bjorken sum rule of the 
free nucleon. However, as in the unpolarized case, the absence of the free neutron target 
means that nuclei have to be used for this purpose. In fact, to study g"^ the SLAC E142 
Collaboration has measured the spin structure function of ^He, while the Spin Muon 
Collaboration (SMC) [|22| has observed the gi structure function of the deuteron. (For recent 
experiments, see Ref. PB[].) 

The nuclear binding, Fermi motion and shadowing modify the parton distributions in a 
nuclear medium. As in the case of Gottfried sum rule in nuclei 0,0, it is worthwhile to 
study how the flavor-nonsinglet combination of the spin structure functions, which appears 
in the Bjorken sum rule, is modified in a nucleus. We can discuss it again using a pair of 
polarized mirror nuclei - the promising pairs may be (^He, ^H) and (^Li, ^Be) (we 
will discuss several candidates in Sec. If the spin structure functions of unstable nuclei 
were measured with a certain accuracy, one could obtain significant information on the spin 
structure of partons in nuclei. 

In general, the spin-average DIS cross section, a, and polarization asymmetry. Act, of a 



target nucleus with mass Mt and spin J in the scaling limit is given by p6 



dxdy 



Cx[l + {l-yf] Y: iFi{x)pi, (1) 

L=even 

-Cx[l~{l~yf] Y: i9i{x)pl (2) 



L=odd 



with a kinematical factor C = c'^MtE /2TiQ'^, where E is the incident lepton energy, q^{= 
—Q^) the four momentum transfer squared and y the ratio of the energy transfer to E in the 
target rest frame. '[Ei and "[gi are, respectively, multipole spin-average and multipole spin- 
dependent structure functions, which are defined in terms of the target hehcity-dependent 
structure functions, E^^ and g(^ {H stands for helicity), and a projection matrix, M/: 

iFi=tr[M/F/], ig, = tT[Mig(]. (3) 

Here Ef and g( are {2J+ 1) x {2J+ 1) diagonal matrices with their elements {EOhh = E(^ 
and {g()HH = gf^ , respectively, and M/ is defined as 

{MDhh' = {-y-\JH, J - H\m)6HH', (4) 

where {JH, J—H\LO) is the Clebsch-Gordan coefficient. In Eqs. (|T]) and (|^) p'l is a multipole 
projection of a spin density matrix, p"^. The multipole structure functions are irreducible 
under rotation. (Note that in the scaling limit '[E2 is related to "[Ei through the Callan-Gross 
relation.) 

For a target with spin J = 1 the DIS cross section is given in terms of three non-vanishing 
multipole structure functions, gFi, \gi and 2-^1, in the scaling limit. This case corresponds, 
for example, to the DIS on ^Li (see Table |I|). The first one is then related to the usual 
spin-average structure functions, Ei, through 



o'^i(^) = ^^t=TE[^r(^) + ^r(^)] = 2V2jtte,{x), (5) 
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where q^^ is the distribution function of quarks with spin s in the target. (Note that all 
flavor labels and charges are suppressed.) The usual F2 structure function can be obtained 
using the Callan-Gross relation. The second one corresponds to the spin structure function 
for a spin-1/2 target, (71, and they are related as \gi = V^gi. The last one is the quadrupole 
structure function for a spin-1 target, which is sometimes called bi [p7| . 

The ^Li target has spin-2 (see Table |). In such a case there are five non-trivial structure 
functions in the scaling limit: qFi, Igi, ^Fi, ^gi and 4F1. Then, qFi, fgi and are analogous 
to Fi, gi and bi, respectively, for J = 1. The structure function '^gi is analogous to the one 
that first appears at J = 3/2, ^^'^gi (see below). The last one, first arises for J = 2. 

The spin of Li isotopes with odd A is J = 3/2 (see Table 1). In this case there are four 
non-vanishing multipole structure functions in the scaling limit: ^o^Fi, '^{^gi and ^i^gi. 

The first one is again proportional to the usual spin-average structure functions, and the 
second one is analogous to the quadrupole structure function, bi, for a spin-1 target. The 
last two are sensitive to quark spin asymmetries, '^{^gi is the analog to gi for a spin-1/2 
target: 

fg^ix) = ^^[^1^%) - 3gr^/^(x) + qf'^^x) - qf'^^x)], (6) 



Q /2 

while 3 (7i is a new asymmetry: 



It is important to keep in mind that all nucleons contribute to the spin-average structure 
functions but that only unpaired valence nucleons mainly contribute to the other multipole 
structure functions. 

In the present paper we concentrate on the spin-average structure functions of ^~^^Li 
and the spin-dependent structure functions of '''^'^^Li, and study those structure functions 
in detail. Measurements of those unpolarized and polarized structure functions could be 
performed by the MUSES project proposed at RIKEN |^ and/or other nuclear facilities in 
the future [|1^,^. The outline of the present paper is as follows. In Sec. |I|, we summarize 
the main features of the lithium isotopes measured in experiments, and calculate their 
wave functions using Hartree approximation. The spin-average structure functions of Li 



are presented in Sec. III. Moreover, we discuss how the neutron structure function can 



be extracted from Li nuclei. The present way will provide a new alternative to study the 



neutron structure function using radioactive ions rich in neutrons. In Sec. |V| we calculate 
the spin structure functions of the ^'^'^^Li isotopes. The multipole spin structure functions 
of ^i^gi and ^3^(71 for a spin-3/2 target are discussed in detail. The effect of the nuclear 
binding and Fermi motion on ^{^gi is less than 10% for x < 0.7, but it becomes important at 
large x, as in the EMC effect. The spin structure function of ^{^gi is novel and it first arises 
for a spin-3/2 target. It is negative for small x but it turns to be positive for 0.2 < x < 0.5. 
For large x it is again negative and the magnitude becomes large because of the Fermi 
motion. We summarize our results and discuss several candidates for testing the Gottfried 
and Bjorken sum rules in mirror nuclei in the last section. 
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II. STRUCTURE OF LI ISOTOPES 



We show the properties of the ®~^^Li isotopes in Table |. It is quite remarkable that the 
two- neutron separation energy, S2n, decreases dramatically as the atomic mass number, A 
= Z + N, increases and that at the neutron drip line S'2„ of ^^Li reaches about 300 keV. The 
one-neutron separation energy, Sn, also diminishes with increasing A (although it shows a 
dependence on the shell structure), while the one-proton separation energy, Sp, increases, 
mainly due to the attractive p-n interaction in nuclei. It is also notable that the rms radius, 
r^, of ^^Li is about 3.1 fm, which is significantly larger than A^^^ ~ 2.2 fm. Furthermore, 
we list the energy levels of the first 1/2^ states of ^Li, ''Be, ^Li, and ^Be in Table |I| [^0, 
as these are expected to provide information on the spin-orbit force for the nucleons in the 
p-shell. 



A. Model for Li isotopes 

During the last 15 years the structure of light /5-unstable nuclei has been studied inten- 
sively. In particular, ^^Li, which is located at the neutron drip line, has attracted a lot of 
attention. There two valence neutrons are so weakly bound to ^Li that they form a halo 
structure. From this point of view, some simple cluster models such as the di-neutron clus- 
ter model [^, where the valence neutrons are correlated strongly enough to be regarded 



as one particle, were first proposed, and those simple models succeeded in reproducing the 
qualitative features of ^^Li. With increasing observed data on ^^Li, it was, however, realized 
that these simple cluster models could not explain the measured experimental data. Thus, 
some authors started calculating the properties of ^^Li using a three-body model which 
was supported by the fact that ^^Li is a Borromean system (^^Li = ^Li +n + n), where all 
the two-body sub-systems are unbound. At present the three-body model is one of the best 
ways to describe the observed properties of ^^Li. However, speaking strictly, the three-body 
model is only valid for structureless particles. In this sense a standard model to describe 
^^Li has not yet been established. 

Our aim in the present paper is to investigate the nuclear structure functions of the Li 
isotopes and, in particular, to determine the momentum distribution of each nucleon in these 
nuclei. Thus, the three-body model is not suitable for our purpose. In addition to this, since 
we consider that the separation energy is more significant for calculating nuclear structure 
functions than the other properties of a nucleus, in the present paper we use a simple shell 
model with Hartree approximation to calculate single-particle wave functions of Li isotopes. 
Note that it implies that the valence neutrons correlate rather weakly in the halo of ^^Li. In 
subsection |1I C| , we will discuss effects of residual two-body interactions. 

In Hartree approximation the total Hamiltonian of a A-nucleon system is given in terms 
of a sum of nucleon kinetic energy and an effective one-body potential, U: 



A 



1=1 



h 



2 



2M 



(8) 



where M (= 939 MeV) is the nucleon mass, and denotes position r^, spin a,, and isospin 
fi of the i-th. nucleon. The A-body wave function, \E'(^i, ■ ■ ■ , ^j, ■ ■ ■ ^^), is then given by 
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the product of single-particle wave functions, (f>ai{Ci) («« stands for all quantum numbers 
associated with the i-th nucleon). In this model the single-particle separation energy is 
simply given by — e^. 

Determining the consistent isoscalar and isovector mean field potentials is very important 
for discussing not only neutron-rich nuclei but also neutron stars [^. However, it would 
need several elaborate experimantal and theoretical works to fix both the potentials uniquely. 
In the present paper, therefore, we suppose as a simple and effective potential as possible. 
It is spherical and consisits of the nuclear, f/^r, and Coulomb, Uc, parts: 



U{A, Z, r,; r) = Un{A, Z, r,; r) + Uc{A, Z; r) 



1-r. 



(9) 



where is the 2;-component of isospin (r^ = +1 for neutron (n) and —1 for proton (p)). 
Here the nuclear part is given by [ P^ ] 



UN{A,Z,Tz;r) 



1 + expf^j '^^'^ 



-f/ 



LS 



1 + exp ( ^ 



(10) 



with R = tq ■ {A — ly/^. In Eq. (|ToD the first and second terms are the central and spin- 
orbit forces, respectively, and we further suppose that the potential depths, Uq and Ulsi 
depend on A, Z, and r^. The geometrical parameters, tq and a, are taken to be common to 
both the protons and neutrons, in order to reduce the number of parameters for the single- 
particle potential. The last assumption is also based on the fact that the ^~^Li isotopes are 
considered to have the normal density distribution (see Table 2 and Fig. 1 of Ref. f^). (As 
will be shown in Fig. ||, however, the present potentials yield decoupled proton- and neutron- 
density distributions in ^Li and ^Li nuclei. We will discuss the reason for this discrepancy 
in subsection |11 C] .) 

Concerning the Coulomb part, we assume that it is formed by a Gaussian charge distri- 
bution with the charge oi Z — 1 and a radius parameter Oc, 



Pe(r) = (Z-l)e- (aev^)-3e 



(11) 



We then obtain 



Uc{r) 



[Z-l)e^ 



(12) 



where erf (a;) is the error function and Oc is determined from the rms radius of the 
Gaussian distribution 



(13) 



We here use the experimental rms radii of the matter distributions of the Li nuclei except 
for ^^Li. For ^^Li, we use the matter distribution of ^Li as the rms radius of the charge 
distribution. 
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B. Shell configuration 



Assuming that identical nucleons in a given shell orbit are coupled to a spin-singlet state 
and that no two-body residual forces exist, we can fix the configurations for the ground 
states in the Li isotopes, based on the simple shell model. For example, one can write the 
shell configuration for ^Li as 

(7rlSi/2)^(7rlp3/2)^(z/lSi/2)^(z/lp3/2)^ • 

We also assume that the valence neutrons in ^^Li occupy the i'lpi/2 orbit although higher 
configurations such as and 2si/2 are thought to be involved in forming a realistic halo 
structure. In our model the single-particle excitation energies of a neutron and proton are 
given by 

Ae{n) = I enilpi/2) - Sn{lp3/2) I, A£:(p) = | ep{lpi/2) - ep{lp3/2)\, (14) 
respectively, and the rms radii of ^~^Li and ^^Li nuclei are given by 

^(2(r').i.,/, + 2(r2).i.,/, + {r').^^.,,, + {A - 5)(r2).i,3/, ), (15) 

=11 = ji^2{r')^u,/, + 2{r'),u,/, + (r').ip3/. + 4(r^).ip3/. + 2(r2).ip,/, ) , (16) 



r-2 

' A=6-9 



respectively. 



C. Model parameters and comments 



In order to obtain the potential parameters, we divide the Li isotopes into two groups: 
one consists of ^Li and ^Li, and the other consists of ^Li, ^Li and ^^Li. In the former and 
latter groups we start with ''Li and ^Li, respectively. The energy differences between the 
first 1/2^ (lPi/2) and ground states (1^3/2) in ^Li and ^Be are less than 500 keV, while 
those in ^Li and ^Be are less than 3 MeV (see Table |T[). We therefore assume that the 
spin-orbit force in '''^Li is not so strong and that it does not affect much the rms radius 
and single-particle energy levels. First, for each of ^'^Li we roughly determine four potential 
parameters, r^, a and Uq, for n and p, so as to produce the rms radius, r^, as well as the 
one-neutron and one-proton separation energies, Sn and Sp, with a constraint that the four 
parameters obtained are not far from the standard values of the potential parameter set 
given in Ref. |Q. (Note that the two geometrical parameters, tq and a, are common to 
both p and n.) 

We next added the spin-orbit potentials, Uls, for ^ and p to fit approximately the energy 
difference, Ae, between the first excited and ground states. Then, we finally tune all the 
potential parameters to exactly reproduce Sn, Sp, Ae{p), Ae{n) and in the vicinity of 
the values determined in the rough fitting. For the other nuclei in each group, we start 
parameter-fitting using the values for ^Li or ^Li and find an appropriate set under the 
condition that the parameter values differ from those for ^Li or ^Li as little as possible. 



We present the optimum parameter set in Table [II and show how the present model 



reproduces the experimental data in Table IV. We have also checked numerically how the 
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properties of the Li isotopes depend on the choice of parameter values. It is then found that 
the nuclear properties depend on Tq and Uq somewhat strongly, but that its dependence on 
a and U^s are weak as we expected first. 

We comment here on the one-neutron separation energy of ^^Li in Table |ITT[ The value 
of 5*^ obtained from the experimental binding energy is 0.326 ± 0.042 MeV, which is larger 
than its two-neutron separation energy, 0.301 ± 0.029 MeV (see Table |). Therefore, the 
probability for two-neutron removal may be larger than that for removing one-neutron. 
(Note that ^°Li is observed experimentally as a resonance state.) Thus, in the present 
calculation we suppose two possibilities for the one-neutron separation energy of ^^Li in 
order to see the dependence of the structure functions on Sn- We take Sn = 301 or 350 keV, 
and hereafter we call the case of Sn = 301 (350) keV ^^Li-a(b), respectively. 

As the neutron number increases the central potential for neutrons becomes shallow, 
while that for protons becomes deep (although they show an odd-even mass number effect). 
This is vital to fit the small separation energy of the valence neutrons as well as the observed 
rms radii of Li isotopes. Remarkably, the strength of the spin-orbit force for both protons 
and neutrons grows drastically with increasing A. In Table [V| we present the energy levels 
of the single-particle states in the Li isotopes. The rms radius calculated for each single- 
particle wave function is also given. According to the weak binding of the valence neutrons 
in Li, the rms radius of the last neutron becomes large with growing A. 

Here we should point out some caveats for the present model (see Table |): (i) In the 
observed data the two- neutron separation energy, S2n, is not twice the one- neutron sepa- 
ration energy, while the present calculation gives S2n = '^Sn- (ii) In the low- lying spectra 
of ^Li and ^Li their degeneracies are decoupled in the experiments, (iii) The distributions 
of protons and neutrons in ®'^Li are quite different at large r in the calculation (see also 
subsection [11 D| ). Those difficulties indicate the necessity of including residual two-body in- 
teractions in nuclei. The first defect could be amended by considering the effect of a pairing 
force between neutrons, while the last one may indicate that a p-n pairing is required to 
suppress the leakage of neutrons in ^'^Li. Note that neutron leakage or shrunkage of the 
proton distribution infiuences the nuclear structure functions much less than the separation 
energy. This fact allows us to use the present parameter sets for their single-particle poten- 
tials, where tq and a are common to both p and n, in spite of the caveat (iii). Furthermore, 
since we used effective potentials which reproduce the separation energies of nucleons close 
to the Fermi surface, and it is hard to measure experimentally the separation energies of 
deeply bound nucleons, it may not be necessary to be too concerned about the caveats (i) 
and (ii) in the present calculation. 



D. Nucleon kinetic energy and average binding energy 

In Figs. |I] - H we present the wave functions in coordinate space, MQ,(r), and the radial 
momentum distribution of a nucleon, (paik), in the Li isotopes. One can see that as the 
neutron number increases the tail of the neutron wave function in the Ip orbit is longer. 
Accordingly, the distribution of the Ip neutron becomes sharp in momentum space. As 
expected, the valence neutrons in the halo of ^^Li distribute very broadly in coordinate 
space and their momentum distributions are quite sharp around /c ~ 0.1 fm~^. The density 
distributions are illustrated in Fig. ^ As pointed out in subsection |11 (J| , the neutrons leak 
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out to large r even in ^'^Li, which may imply that a proper (attractive) p-n pairing force is 
required. 

In an unpolarized nucleus the kinetic energy of a nucleon in the single-particle shell, 
a = {n,£,j,Tz}, can be calculated by 

it^) = ^iT. j ^,\<Pn,.rnrAk)? ^ = ^ (17) 

Here (pnejmr^ik) is the nucleon momentum distribution, which is given by Fourier transfor- 
mation of the wave function in coordinate space, and it is normalized as 

I £k\(l)nejmrAk)\^ = 1. (18) 



(27r)3 

The kinetic energy of the center-of-mass (cm) motion of the nucleus is then evaluated by 

Tcm = jY:{taJ, (19) 

and the spurious cm motion should be subtracted from the total energy. We thus obtain 
the mean- nucleon kinetic energy, (t), in the nucleus and the total binding energy: 

(^) = ^(i-^)E(U, (20) 

(5) = ^E(^^ + (^-.))-^-, (21) 

respectively. In Table |VT| we summarize the kinetic energy of the nucleon in each shell model 
orbit. We should note that the calculated binding energies are rather smaller than the ob- 
served ones (see Table |). With increasing the neutron number the total binding energy 
decreases. Probably this discrepancy may be improved by taking residual two-body corre- 
lations into account. However, we again emphasize that the present calculation reproduces 
the observed separation energies of Li isotopes, which plays the most important role in the 
calculation of nuclear structure functions. 



III. SPIN-AVERAGE STRUCTURE FUNCTIONS OF LI ISOTOPES 

A. Light-cone momentum distributions and nuclear structure functions 

In the DIS of high energy charged leptons off an unpolarized nucleus one can observe 
the distribution of the longitudinal momentum fraction, x, carried by the quarks in the 
nucleus [|]. The nuclear structure function, F^{x), is usually described by a convolution of 
momentum distributions of nucleons in A and the free nucleon structure function pl,^: 

F^'^ix) = E rdyUMy)Fi{x/y), (22) 
i -^^ 
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where F2 = '■"'^ when the i-th nucleoli is a proton (neutron). (Note that in the scahng hmit 
F^{x) is related to F{^(x) through the Callan-Gross relation. See also Eq. (^.) fajAiv) is a 
distribution of the longitudinal momentum fraction, y, carried by the nucleon in the nucleus 
A, which is defined as 

0<2/M=-^<l {0<y<A), (23) 

with P'^ the four momentum of the nucleus A. 

In the scaling limit the light-cone momentum distribution in the target rest frame is then 
written as 

/.Mfe)-/(f,S«,MW(l4)<.-^^). (24) 
with the nucleon spectral function 

S^jAik) = 27r6{ko -M-e, + Tn)\^aXk)\\ (25) 

where — £«(> 0) is the separation energy of the z-th nucleon in the present approximation, 
Tr = k"^ /2Ma^i{> 0) the recoil energy of the residual nucleus {A — 1) with its mass Ma-i, 
and "^a-ik) the three-momentum distribution of the struck nucleon i. For an unpolarized 
target, |\E'q,. (fc)p is given in terms of the momentum distribution of the nucleon as 

\^c.Xk)\' ^ 77^T.\<l^nejmrAk)\' = ^\Mk)\^ (26) 

which is spherically symmetric and independent of the angular variables. 
Inserting Eq. ( pS]) into Eq. (^4]), we find 

U/A{y)= J (^2vr5(fco-M-5, + T«)|vl/„^(A;)r(l + -j5(y-^— j. (27) 
Using a new notation 

M = M{k) =M + ei-TR, (28) 

we can rewrite Eq. ( P7| ) as 

d^k ^ k^\/ k° + k^^ 



UMV) = / (^2rf(*„ - A#)|*„.(*)|^(l + ^S{y 



2(27r)3 

The distribution fai/A{y) vanishes unless 



M / M 



< 1, (30) 
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that is, 



\My-M\ = \My-M -ei + TR\<k. (31) 

Since the recoil energy, Tr, depends on A;^, Eq. (^) yields a transcendental equation. 
However, Tr is practically small and hence we may replace Tr by its expectation value 

Defining the lower limit of the fc-integral in Eq. (^) as 

kL = \My-M-e, + {TR)^^\, (33) 

we find 

/../.M-^.£^K.(*)l^ (34) 

The longitudinal momentum distribution, fai/A{y), must satisfy the condition of baryon 
number conservation: 



A 

dy fci/A{y) = '^- (35) 



We have confirmed that fat/Aiy) obtained numerically satisfy the condition given in Eq. (|35D 
with an accuracy of better than 0.2 %. In Figs. ^ — ^ we present the f{y) distributions for 
^~^^Li. Remarkably, the halo neutrons in ^^Li are distributed quite sharply and the shape 
of f{y) is almost symmetric around y = I. This means that to a good approximation such 
neutrons can be regarded as being essentially free. 



B. Numerical results for 

Here we present numerical results for the spin-average structure functions of the Li 
isotopes, which were calculated using Eqs. P^ ) and (|3^). For the structure function of 
the free nucleon, the CTEQ5 parameterization |^ was used. Fig. § illustrates the EMC 
ratios of the Li isotopes, which are defined by {F^{x)/A)/{F2{x)/2). We should note that 
the nuclear structure functions in the EMC ratio are usually corrected to be iso-scalar - 
the so-called iso-scalarity correction However, we did not apply the correction to the 

present ratio, and hence the Li structure functions in the numerator are not iso-scalar except 



for Li. The structure function of the deuteron was calculated using the Paris potential [|17 . 
In order to evaluate its recoil energy, (Tr), we employed relativistic kinematics, since the 
residual debris is just one nucleon. We find that (Tr) ~ 9.0 MeV and —e + (Tr) ~ 11 
MeV [^. (If one chooses nonrelativistic kinematics (Tr) is about 9.4 MeV ||36|| .) 

In Fig. H each EMC ratio increases rapidly with x in the region x > 0.7 and has a 
minimum at around x = 0.6. Also the minimum grows deeper as the number of neutrons 
increases. The rapid increase of the EMC ratio is caused by the Fermi motion, which 
can enable the nuclear structure function to have a tiny but finite component at a; > 1. 
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Concerning the depression from unity appearing at intermediate x, one needs a careful 
discussion because no iso-scalarity correction is applied in the present EMC ratio. 

For ^Li, the depression from unity around x = 0.6 is rather small. From Tables |V] and|V| 
one can calculate the average-separation energies of protons, neutrons and all the nucleons, 
as well as the average- recoil energy for each Li isotope. They are presented in Table |V11 . 
The average proton- and neutron-separation energies, {e)p and {e)n, of ^Li are about 15 
MeV and 17 MeV, respectively, while the average-recoil energy (Tr) is about 2.3 MeV. The 
nuclear binding effect, which is given by the sum of them, is thus about 18 MeV for ^Li, 
while it is about 11 MeV for the deuteron. The small difference between these values is the 
reason for the small depression in the EMC ratio of ^Li. 

For /3-unstable Li isotopes, which are not iso-scalar, the minima appearing around x = 0.6 
in the EMC ratios become deeper as the neutron number increases. Let us consider the case 
of ^^Li as an example. From Table ^ we find = 6.0 (5.9) MeV, {e)p = 30 (28) MeV, 



the average-separation energy of the total system, {e)tot = \{N{e)n + Z{e)p) = 12 (12) MeV, 
and (Tr) = 1.4 (1.3) MeV for ^^Li-a(b). Thus, the nuclear binding of ^^Li is about 13 MeV, 
which is much closer to that of the deuteron, compared with the case of ^Li. It is then 
expected that the depression in the EMC ratio of ^^Li would be more suppressed than that 
of ^Li. However, the EMC ratio shown in Fig. ^ is completely opposite. Since /3-unstable 
Li isotopes have more neutrons than protons, as the neutron number increases the neutron 
structure function should contribute more to the EMC ratio than does the proton. 

Let us roughly estimate the effect of neutron excess on the EMC ratio. We suppose 
that a nucleus is a simple collection of protons and neutrons and that the binding effect is 
ignored. It is well known that -Ffl^) > F^ix) and around x = 0.6 the structure function 
of the free neutron is about a half of that of the free proton [^]. In this case the EMC 
ratio of ®Li to the deuteron is just unity because ^Li is iso-scalar. Using the relationship 
Fi{x) ~ 2F^{x) one finds that the ratios for 7,8,9,iiLj become 0.95, 0.92, 0.89 and 0.85 at 
X ~ 0.6, respectively. It is remarkable that those values are already close to the calculated 
ratios at x = 0.6 in the figure. When the effects of the nuclear binding and Fermi motion are 
turned on, they, of course, suppress the Li structure function at intermediate x and enhance 
it at large x. However, such a suppression at intermediate x may not be large because it is 
approximately determined by the difference between the sum of the average separation and 
recoil energies of Li and that of the deuteron - see Table |VI1| . 

If those successive EMC ratios could be measured in the future one would obtain system- 
atic information on the off-mass shell effect on the neutron structure function in a nuclear 
medium flM- We also note that the EMC ratio of ^^Li is not sensitive to the choice of the 



one-neutron separation energy, S'„ (see Fig. |^). 

We next discuss how the last neutron in each Li isotope is bound as compared with the 
neutron in the deuteron using the following ratio [0 

F2(^Li) -F2(^-iLi) , , 

where F2(^Li) is the nuclear structure function of "^Li. Note that for Rdill) take the 
numerator to be [F2(^^Li) — F2(^Li)]/2. The results are plotted in Fig. In the region 
X < 0.7 the ratios Rd{A) exceed unity, except for A = 7. This behavior may reflect the 
tendency for the energy level of the neutron- lp3/2 orbit to decrease with increasing neutron 
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number (see Table 0), even though the odd-even mass number effect can be seen in those 
energy levels. In large x region Rd{A) becomes less than unity, which is consistent with the 
fact that the momentum distribution of the last neutron in the Li isotope does not have a 
large high- momentum component. However, one has to keep in mind that the contribution of 
the nucleons in the core nucleus to the structure function F2('^Li) is not completely cancelled 
by F2('^~^Li) and that there remains an unavoidable contamination of the binding effect on 
the proton in the deuteron in the denominator. 

Let us turn to the discussion of the following ratio 

^,.) . ^-(-") (37) 

It would be useful to extract the neutron structure function from data for the series of Li 
isotopes. (For A = 11 we set the same numerator as in Rd{ll).) If one uses the ratio, 
Rn{A), and the observed structure functions of ^Li and ^~^Li, then one can obtain the 
neutron structure function from [F2(^Li) — F2{^~^Li)]obs./ Rn{A) . Figs. |10| and |TT] show 



Rn{A) for Li isotopes as well as that for the deuteron, Rn{D) = (F^ — F^/F^. In the case 
that the neutron structure function is assumed to be simply given by F^ — Ff', one can see 
from Rn{D) that the correction due to the nuclear binding is about 4 % at around x = 0.5 
and that the effect of Fermi motion is quite large in the region x > 0.7 (about 10 % at 
X = 0.7). 

On the other hand, some combinations of the Li structure functions provide a pretty flat 
ratio below x ~ 0.6 — 0.7. Within the present model, for example, the ratio -Rn(9) is almost 



unity for x < 0.5 and exceeds unity by only about 2 % at x = 0.6. We find in Table |VI1 
that the average-separation energies as well as the average-recoil energies of ^Li are almost 
identical to those of ^Li. This is the reason why -R„(9) is so close to unity at small and 
intermediate x. The ratios -Rn(ll — and -Rn(ll ~ b) are also fiat below x < 0.5 — 0.6 
because of the weak binding of the neutrons in the halo of ^^LiQ (See also Table |ITT| .) The 
ratios are about 1.03 at x = 0.5. In particular, in -R„(ll — b) the deviation from unity, 
\Rn{ll -b)- 1.0|, is below 3 % for x < 0.75. 

However, in the large x region the ratio increases in -Rn(ll — a), while it decreases in 
i?„(ll — b). For the case Sn = 301 keV (^^Li-a) the valence neutrons in ^^Li are bound more 
loosely than those for the case of Sn = 350 keV (^^Li-b). Correspondingly, the rms radii 
of the neutrons in the lpi/2-orbit are, respectively, 4.91 fm and 4.82 fm for the former and 
latter cases (see Table 0). In order to reproduce the measured rms radius of riiLi = 3.10 
fm, the nucleons in the core need to be more tightly bound for the former case than for the 
latter (see again Table 0). Thus, the effect of Fermi motion due to the nucleons in the core 
appears to be stronger for the former case than that for the latter, which makes -R„(ll — a) 
increase and -R„(ll — b) decrease in the region x > 0.7. From the behavior of -R„(ll) at 
large x, we can say that the ratio is quite sensitive to how loosely the neutrons in the halo 
of ^^Li are bound, and that the ratio strongly depends on high momentum components in 
the nucleon momentum distributions for ^Li and ^^Li. It is certainly intriguing to work 



^Recall that the cases "^^Li(a-b)" correspond to the two choices of separation energy introduced 
earlier, Sn = 301 (a) or Sn = 350 (b) keV. 
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out a more elaborate calculation for the nuclear structure of the Li isotopes, where residual 
interactions such as pairing correlations are taken into account, in order to study the ratio 
Rn{A) in detail. 

We have already noticed that the energy spectra of the Li isotopes show the odd-even 
mass number effect (see Table |V|). As a method for cancelling the effect simply, it would be 
interesting to see the following ratios: 

, Fsf^Li) -F2(^-2Li) , , 

Rn2{A) = > i. (39) 

Fig. ig presents the ratios for ^Li and ^Li. (Note that Rd2{,A = 11) = /?d(ll) and i?„2(ll) = 
i?„(ll).) One can again see that the ratios for a: < 0.6 are pretty close to unity and the rapid 
increase due to the Fermi motion appears for x > 0.6. In particular, the deviation of i?d2(8) 
(i?„2(9)) from unity lies within 2% (3%) for x < 0.55. Interestingly, Rd{9) decreases in the 
large x region, while Rd2{^) increases there. It means that the ^Li nucleus has relatively 
larger high-momentum components than ^Li, which may partly result from the fact that 
the binding energy of protons in the lsi/2 orbit for ^Li is much larger than that for ^Li (see 
Table [Vp. 



IV. SPIN STRUCTURE FUNCTIONS OF LI ISOTOPES 



Contrary to the case of unpolarized structure functions , so far, only a few theoretical 
studies have been made for the effect of nuclear binding on the spin-dependent structure 
functions 



14,17,25,38-42 



The present situation seems to be that the errors in extract- 
ing nuclear effects from the spin- dependent structure functions still seem to be within the 
experimental error bars. However, recent progress in high precision measurements of spin- 
dependent structure functions for the neutron, proton, and deuteron [ P0| , p3[ , together with 
improved parameterizations for the polarized parton distribution functions (PPDF) 
will make it necessary to know the effect with a high accuracy in near future. 



3,44 



A. Nuclear spin structure functions 

Here we systematically consider the spin structure functions of the Li isotopes. Let us 
suppose that the spin-dependent structure functions of a nucleus separate into the contribu- 
tions from those of protons and neutrons. In the conventional nuclear model the polarization 
of a nuclear ground state is determined by the unpaired nucleons because two identical nu- 
cleons in the same shell orbit prefer coupling to spin of zero. Thus, the spin of a Li isotope 
with an odd mass number is 3/2, which is attributed only to the valence proton in the lp3/2 
orbit, while the spin of the ^^^^ Li ground state is 1 (2), which is given by a combination of 
the unpaired proton and neutron in the lp3/2 orbit (see Table |). 

It is, however, important to consider residual interactions between single-particle states 
in a nucleus, which lead the configuration mixing and hence a change of the wave function of 
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the valence nucleon. If we ignore the nuclear binding and Fermi motion effects the nuclear 
spin structure function may be approximately given by |25,40 



gt{x) = X + X gl{x), (40) 

where Pp{n) is the effective polarization of proton (neutron) in the nucleus A. It is well known 
that, for example, in ^He the tensor force etc. provide a mixture of S'- and D-states to the 
S-wave state and hence that the neutron is somewhat depolarized and the proton gains a 
slight downward (negative) polarization due to the configuration mixing: P^^^ = 86 ± 2% 
and Pp""^ = 2 X (—2.8 ± 0.4)% ^0|. In the series of Li isotopes a similar situation would 
occur because of a mixing of higher configurations. It is therefore necessary to perform a 
calculation of the nuclear structure of Li isotope including residual interactions to get precise 
values of the effective polarizations. We do not perform such a calculation in this paper. 

Instead, in this section we consider the effects of nuclear binding and Fermi motion on 
the spin asymmetry for a target with spin 3/2 in detail. This corresponds to DIS off "^'^'^^Li 
targets. For such targets there are four non-vanishing multipole structure functions in the 
scaling limit: ^o^-Fi, ^2^Pi, ^^iQi and ^f^Qi - recall the discussion in Sec. |. The first one is 
proportional to the usual spin-average structure functions, which has been studied in the 
previous section, while the second one is analogous to the quadrupole structure functions, 6i, 
for a spin-1 target. The last two provide information on quark spin asymmetries: { gi is the 
analog to gi for a spin-1/2 target, while ^i'^gi is a new asymmetry function. It is important 
to keep in mind that in the present approximation only unpaired valence nucleons contribute 
to the multipole spin structure functions. 

B. Spin structure function for a spin-3/2 target 

The spin asymmetry for a spin-3/2 Li target is given by 

= Cx[l-[l~yr\ 

L=l,3 



^^^y Cx[l-{l-yr] J: %'g,ix)pr, (41) 



where the factor C is given below Eq. (^. Then, we can calculate the multipole spin 
structure function using the following convolution form: 

%i(x)= rdy-gf{y)g{{x/y), (42) 
Jx y 

where g^^'^iy) is the light-cone spin distribution given by tr[Mf/'^ gi^'^] - note that gi^'^ is a 
diagonal (2 J + 1) x (2 J + 1) matrix defined below Eq. (H). 



We find g^"^ exphcitly as 



9l^\y) = ^co(i/) - ^C2(y), (43) 

9l%) = -^^^(y)^ (44) 
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where 



2'K6{ko 



X Pl(cos9)5 y 



M j ■ 



A;0(A;0 + M) 



(45) 



In Eq. (|45|) the term of the square blacket is the so-called flux factor, which stems from 
the matrix element of 7''' 75 with respect to the nucleon spinor. In the limit — M it 
coinsides with the flux factor given in Refs. P,^. Instead of the present flux factor, the 



factor, (1 + k^/k^), which is the same as that for the f{y) distribution, is often used to 
calculate the nuclear spin structure function. However, the difference between two cases is 
less than 2% in the y distribution (see Fig. [T^ in subsection |1V (Jl ). 

After performing the integrals with respect to k^ and the angular variables we finally 
obtain 



co{y) 
C2{y) 



kdk 



M2 

r'^ kdk 



M JkL 4(27r)= 



fc2 



2M2 



,M2 



y- 



M 



- 1 



y + ^^y 



1?- 



2M2 



(46) 
(47) 



where <Pip^^2{k) is the radial wave function of the valence proton in momentum space. The 
term of 0{e/M) is neglected in the flux factor because its contribution is very small. 



C. Numerical results for spin structure functions 

In Fig. |l3|we plot the distributions, co{y) and C2{y), for ^'^'^^Li isotopes. For comparison. 



we also present the light-cone momentum distribution, f{y), of the valence proton in the 
lps/2 orbit. One sees that the shape of the Co{y) distribution is very similar to the f{y) 
distribution, while that of C2{y) is quite different. Note that the small difference between 
Co{y) and f{y) is caused by the second and third terms of the flux factor in Eq. (|^). It 
should be noted that CL{y) for ^^Li is not sensitive to the choice of the one-neutron separation 
energy, Sn, and hence that the difference between the ciiy) distributions for Sn = 301 and 
350 keV is quite small. 



There is a sum rule for the multipole spin structure function ^6 

rA 
10 



f dxx"^-^ igi{x) =0, for 1 < n < L. (48) 
Jo 

This offers a sum rule for C2{y) 

[^dyc2{y)=0, (49) 
Jo 

while Co{y) is free from such a constraint. We have checked the sum rule for ^s^gi{n = 1) 
numerically, and found that the absolute value of the integral, Eq. (ff3), is less than 0.0014 
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for the ^'^'^^Li isotopes, where the calculation was performed in the region 10~^ < a; < 2.20. 
The sum rule is certainly satisfied within the numerical precision. 

In Figs. and |1^ we present the multipole spin structure functions, ^^iQi and ^3^(71, 
for A = 7, 9, 11 normalized Qi for the free proton. For the numerical calculation we have 
used the polarized parton distribution functions (PPDF) parameterized by D. de Florian 
et al. g and M. Gliick et al. H (denoted by GRSV PPDF)_. We use one of the options 
for both, up to the next-to-leading-order (NLO) within the MS scheme, which includes the 
effects of flavor asymmetries in the proton sea. The calculations are performed at = 10 
GeV^. Since the two parameterizations eventually gave very similar results we show only 
the result with the GRSV parameterization in the figures. Moreover, since the choice of the 
one-neutron separation energy, 5'„, in ^^Li does not affect the calculated result much, we 
illustrate only the ratios for ^^Li-a (i.e., Sn = 301 keV). 

In graphing the ratio for '^{^gi (see Fig. 0) we have divided out the "geometrical factor" 
2v^/3 so the plotted ratio would be unity in the absence of the nuclear binding and Fermi 
motion. The shape of the ratio is very similar to the usual EMC ratio - see also the 
distributions of f{y) and co{y) in Fig. 13. The effect of the nuclear binding and Fermi motion 
on '^{^gi is then about 10% for x < 0.7, while the Fermi motion effect is very important for 
large x. This observation is consistent with the earlier results for the deuteron |]4l|j4^ and 
^He |E0|. For x < 0.6 the ratio for ^Li is a little shallower than that for ^'^^Li. This is 



because the weak binding of the valence proton in ^Li as compared with in ^'^^Li isotopes - 
see Table 0. 

The shape of the ratio for ^^^gi is quite different from that for ^{"^gi- It is negative for 
small X but it becomes positive for 0.2 < x < 0.5. However, the magnitude is very small in 
the region of x < 0.6. For very large x it is again negative and its absolute value becomes 
large because of the Fermi motion. Because of the weaker binding of the valence proton in 
^Li than in ^'^^Li, the magnitude of the ratio for ^Li for x > 0.5 is smaller than that for 
^'^^Li. We should note here that the present result is quite different from that of Jaffe and 
Manohar [|26| . In particular, the ratio for large x is negative in our realistic calculation, while 



it is large and positive in Ref. This discrepancy is caused by an approximation used 
in Ref. [26|, where the distributions of Co(?/) and C2{y) are expanded in terms of 5-functions. 
This amounts to neglecting the positive wings of 02(2/) in Fig. which is numerically 
inaccurate. 



V. SUMMARY AND CONCLUSION 

We have studied the spin-average and spin-dependent structure functions of the lithium 
isotopes. First, we have calculated the wave functions within Hartree approximation. In 
the present calculation we used a phenomenological nuclear potential of the Woods-Saxon 
type, in which parameters are fitted so as to reproduce the separation energies of protons 
and neutrons, the single-particle excitation energies and the rms radius for each Li isotope. 
In the calculation, for ^^Li the separation energy for the last neutron in the halo was taken 
to be either Sn = 301 or 350 keV in order to examine the S'„-dependence of the structure 
function. Next, using the wave functions we have individually calculated the light-cone 
momentum distribution of the nucleon in each shell model orbit. We have found that the 
f{y) distribution of the valence neutron in ^^Li is very sharp and almost symmetric around 
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y = 1, which means that this neutron is bound under a condition quite close to the neutron 
in free space. 

Adopting the usual convolution technique we have calculated the spin-average struc- 
ture functions of Li isotopes. The EMC ratios of the Li nuclei were first illustrated and 
the role of the neutron structure function in the EMC ratio was discussed. As the neutron 
number increases the neutron structure function plays a more important role in the nuclear 
structure function because its weight becomes larger than that of the protons. Furthermore, 
we have presented the ratio, R^, which is defined by the ratio of the difference between the 
Li structure functions of mass number A and A — 1 to the difference between the structure 
functions of the deuteron and the free proton. The ratio, -R„, which is defined as in i?^ but 
the denominator is the free neutron structure function, was also illustrated. The ratios show 
the sensitivity of the nuclear structure function to the last valence neutron. In particular, 
the ratios -Rn(9) and are very useful for extracting the free neutron structure function 

from Li data. We have also investigated the ratios Rd2 and Rn2, where the numerators are 
chosen to be F2("^Li) — F2('^~^Li), to eliminate the odd-even mass number effect in the Li 
structure. 

We have studied the spin-dependent structure functions of Li isotopes. In the present 
study we have concentrated on Li nuclei with odd A and calculated the multipole spin 
structure functions, ^{^gi and ^i^gi, where the former is the analog to gi for a spin-1/2 
target and the latter is a new spin asymmetry which first arises for a spin-3/2 target. We 
have studied the nuclear binding and Fermi motion effects on those spin structure functions 
in detail. Here the spin structure function is assumed to be given by a convolution form of 
the light-cone spin distribution of the valence proton, which is provided by a combination 
of Co{y) and C2{y), and the spin structure function of the free proton. 

We have found that the shape of the ratio of ^{^gi to g^ is very similar to the usual EMC 
ratio, while the ratio for 3 gi is quite different from the EMC one. It has been confirmed 
that the effect of the nuclear binding and Fermi motion on ^{^gi is about 10% for x < 0.7, 
but that the Fermi motion effect is very important for large x. On the other hand, the 
multipole spin structure function of '^^'^gi is negative for small x but it turns to be positive 
for 0.2 < a; < 0.5. For very large x it is again negative and the magnitude of the ratio 
becomes large because of the Fermi motion. The present observation contradicts an earlier 
reported result, where an approximation was used to calculate the distributions of co{y) and 
C2{y)- 

We should note that in both the spin-average and spin-dependent cases a shadowing (and 
ant i- shadowing) correction would be necessary for the structure functions at small x{< 0.2). 



This would require delicated calculation along the line of Refs. [p4|,p5|. For further studies 



it will also be necessary to include two-body residual correlations in the calculation of the 
structure of Li. 

Here we would like to comment on the modification of the Gottfried and Bjorken in- 
tegrals in nuclei: the former is defined by an integral of the difference between the spin- 
average structure functions of a pair of mirror nuclei (divided by x) |]7|, pT| , p^ , while the 



latter is given by an integral of the difference between the spin structure functions of mirror 
nuclei p^ , p5| . As mentioned in the introduction, such quantities could provide us significant 
information on the flavor-nonsinglet parton distributions in a nucleus. Together with the 
flavor-singlet structure functions, such information is quite important for determining the 



20 



individual isospin distributions in a nucleus separately. 

We list several candidates for such pairs of mirror nuclei. The most promising candidate 
is of course the pair ^He and '^H [^|TT],|I2|, where both the nuclei are stable under strong 



interactions. However, this is also the lightest pair, so that the nuclear modification on the 
structure functions cannot be expected to be large. The pair ''Li and ''Be should also be 
good because ^Be has a long half-life time (53 days) and ^Li is stable. We thus believe 
that the Li-Be pair may be the best one for testing the nuclear effects on the Gottfried and 
Bjorken integrals in a medium. The next candidates would be the pairs ^^C-^^B and ^^N- 
^^C, where the half-lives of ^^C and ^^N are about 20 and 10 minute, respectively, while ^^B 
and ^^C are stable. For other candidates one can find some pairs, whose half-life times are 
longer than 1 minutes: ^^O (71 sec)-^^C (5730 years), ^^O (120 sec)-^^N (stable), and ^^F 
(65 sec)-^''0 (stable). However, it is nearly impossible to measure their structure functions 
with fixed-target experiments. 

If DIS with high energy electrons (or muons) off unstable nuclei (like the series of Li 
isotopes) is realized by using a collider machine (like MUSES or (Tesla)HERA) in the fu- 
ture P, p^|28| , available mirror pairs would be widely extended and then quark nuclear 
physics could be developed quite systematically. Although we have discussed only the series 
of Li isotopes in this paper, we can point out other interesting radioactive isotopes: for 
example, the series of isotopes of Be, B, C, O, Na and Mg are of great promise 0. In par- 
ticular, the series of oxygen isotopes are interesting. In the series the smallest and largest 
mass (neutron) numbers are 13 (5) and 24 (16), respectively. For neutron-rich isotopes, a 
neutron skin or a neutron halo is expected to be formed around the tightly bound ^^O core. 
If we could control the atomic number and the difference between the proton and neutron 
numbers freely in measuring the nuclear structure functions, it would stimulate a great deal 
of works that should lead to interesting new quark physics in nuclei 0. 
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TABLES 



TABLE L Experimental data for Li isotopes |Qj5|. The one-neutron, Sn{Z,N), one-proton, 
Sp{Z,N), and two-neutron, S2n{Z,N), separation energies of the nucleus {A = N + Z) 
are defined by 5„(Z,iV) = B{Z,N) - B{Z,N - 1), Sp{Z,N) = B{Z,N) - B{Z - l,N) and 
S2n{Z,N) = B{Z,N) — B{Z,N — 2), respectively, where B{Z,N) is the total binding energy of 
the nucleus "^Z. The rms radius of "^Li is denoted by rA- The numbers in the parentheses stand 
for the error bars of the observed data. 









B (MeV) 


Sn (MeV) 


S2n (MeV) 


Sp (MeV) 


TA (fm) 


6Li 


l+;0 




31.995 


5.664 




4.589 


2.32 








(0.0005) 


(0.0505) 




(0.0505) 


(0.03) 


^Li 


3/2-; 


1/2 


39.245 


7.250 


12.91 


9.998 


2.33 








(0.0005) 


(0.0009) 


(0.0505) 


(0.0015) 


(0.03) 


8Li 


2+;l 




41.277 


2.033 


9.283 


12.45 


2.37 








(0.0005) 


(0.0010) 


(0.0010) 


(0.0305) 


(0.03) 


9Li 


3/2-; 


3/2 


45.341 


4.064 


6.096 


13.93 


2.32 








(0.0019) 


(0.0024) 


(0.0024) 


(0.0090) 


(0.02) 


"Li 


3/2-; 


5/2 


45.642 


0.326 


0.301 


15.30 


3.10 








(0.0271) 


(0.0422) 


(0.0290) 


(0.0971) 


(0.22) 



TABLE IL The first 1/2- levels of ^'^Li and ^'^Be nuclei. The values presented are those of the 
first excited states (except for ^Be) and correspond to the single-particle excitations specified in 
the fouth column. The value for ^Be is for the third excited state. The spins, parities, and energy 
levels of the first and second excited states of '^Be are 1/2+ (1.684 MeV) and 5/2- (2.429 MeV), 
respectively. 







Ex (MeV) 






^Li 


1/2- 


0.4776 


AE(7rlpi/2 


- ^lP3/2) 


■^Be 


1/2- 


0.4291 




- i^l?'3/2) 


9Li 


1/2- 


2.691 


A£;(7rlpi/2 


- 7rlp3/2) 


9Be 


1/2- 


2.780 
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TABLE III. Parameters in the effective potentials for Li nuclei. The numbers (those in the 
parentheses) are for neutrons (protons). We suppose that the one-neutron separation energy of 
^^Li is 301 or 350 keV, which is, respectively, denoted by ^^Li-a or b in the potential parameters. 





C/o (MeV) 


r (fm) 


a (fm) 


Uls (MeV) 


Oc (fm) 


6Li 


5L00 


1.38 


0.66 


1.25 






(5L36) 


(1.38) 


(0.66) 


(1.32) 


(1.894) 


^Li 


50.05 


1.38 


0.66 


1.26 






(58.10) 


(1.38) 


(0.66) 


(1.32) 


(1.902) 


8Li 


59.30 


1.00 


0.53 


7.27 






(88.01) 


(1.00) 


(0.53) 


(9.00) 


(1.935) 


9Li 


55.09 


1.05 


0.54 


12.46 






(80.15) 


(1.05) 


(0.54) 


(9.00) 


(1.894) 


i^Li-a 


43.81 


1.08 


0.51 


12.46 






(71.72) 


(1.08) 


(0.51) 


(9.00) 


(1.894) 


"Li-b 


40.38 


1.14 


0.51 


12.46 






(66.60) 


(1.14) 


(0.51) 


(9.00) 


(1.894) 



TABLE IV. Comparison of the calculated results to the experimental data for ^ ^^Li nuclei. 
The numbers in the parentheses are the corresponding experimental values (see also Tables | and 

0)- 





Sn (MeV) 


Sp (MeV) 


Ae{n) (MeV) 


Ae{p) (MeV) 


VA (fm) 


6Li 


5.665 


4.589 


0.427 


0.442 


2.33 




[5.664] 


[4.589] 


[-] 


[-] 


[2.32] 


^Li 


7.246 


10.006 


0.429 


0.476 


2.33 




[7.250] 


[9.998] 


[0.4291] 


[0.4776] 


[2.33] 


«Li 


2.033 


12.449 


1.397 


2.643 


2.37 




[2.033] 


[12.45] 


[-] 


[-] 


[2.37] 


9Li 


4.064 


13.930 


2.778 


2.685 


2.33 




[4.064] 


[13.93] 


[2.780] 


[2.691] 


[2.32] 


i^Li-a 


0.301 


15.305 


2.550 


2.643 


3.10 




[0.326] 


[15.30] 


[-] 


[-] 


[3.10] 


iiLi-b 


0.350 


15.301 


2.747 


2.737 


3.10 




[0.326] 


[15.30] 


[-] 


[-] 


[3.10] 
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TABLE V. Calculated energy levels of single-particle states and the rms radius of the nucleon. 



The numbers arc for neutrons, while those in the parentheses are for protons, vn stands for the 
rms radius calculated by the single-particle wave function of the nucleon. 





1^1/2 




lP3/2 




lPl/2 






-e (McV) 


tn (fm) 


-e (MeV) 


rN (fm) 


-e (MeV) 


rN (fm) 


6Li 


2L92 


1.93 


5.665 


2.91 


5.237 


2.95 




(20.39) 


(1.96) 


(4.589) 


(2.98) 


(4.147) 


(3.02) 


^Li 


23.13 


1.96 


7.246 


2.84 


6.817 


2.87 




(27.11) 


(1.89) 


(10.01) 


(2.65) 


(9.530) 


(2.67) 


^Li 


21.41 


1.73 


2.033 


3.15 


0.637 


4.00 




(38.51) 


(1.51) 


(12.45) 


(2.18) 


(9.805) 


(2.26) 


9Li 


21.94 


1.79 


4.064 


2.84 


1.286 


3.54 




(37.49) 


(1.59) 


(13.93) 


(2.23) 


(11.25) 


(2.29) 


iiLi-a 


17.86 


1.95 


2.851 


3.13 


0.301 


4.91 




(36.55) 


(1.67) 


(15.30) 


(2.27) 


(12.66) 


(2.32) 


iiLi-b 


16.99 


2.02 


3.096 


3.15 


0.350 


4.82 




(34.88) 


(1.74) 


(15.30) 


(2.33) 


(12.56) 


\Z.oo) 


TABLE VL Kinetic 


energy of a 


nucleon in each shell orbit. The integration over 


the wave 


function 


in momentum space was performed in the region 


< A; < 6 fm-^ 








jylSi/2 




I^lPl/2 


7rlSi/2 


7rlp3/2 


^lPl/2 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


6Li 


12.65 


16.59 


16.27 


12.31 


16.08 


15.72 


^Li 


12.26 


17.01 


16.77 


13.11 


19.16 


18.96 


^Li 


15.92 


17.07 


13.73 


20.62 


28.89 


27.32 


9Li 


14.78 


18.46 


14.45 


18.59 


27.23 


26.03 


^^Li-a 


12.50 


15.66 


10.73 


16.73 


25.86 


25.12 


iiLi-b 


11.62 


15.09 


10.43 


15.47 


24.37 


23.76 
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TABLE VII. Calculated average-separation and recoil energies for Li isotopes. {e)p, (e)„ and 
{£)toti= (-^(s)n + Z{e)p)/A) are, respectively, for the average-separation energies of protons, neu- 
trons and the total system. The average-recoil energy is denoted by (Tr)(= {t)/{A — 1)). The 
bottom row is for the deuteron ||3^ . 



{e)p (MeV) {e)n (MeV) {e)tot (MeV) {Tr) (MeV) 



^Li 15.1 16.5 15.8 2.3 

^Li 21.4 15.2 17.9 2.1 

^Li 29.8 9.78 17.3 2.4 

^Li 29.6 10.0 16.6 2.1 

i^Li-a 29.5 5.97 12.4 1.4 

i^Li-b 28.4 5.88 12.0 1.3 

D — — 2.22 9.0 
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FIG. 1. Radial wave functions, Ua{r), and momentum distributions, (j)a{k), of the single-particle 
states in ^Li (left panels) and ^Li (right panels). The top two panels are for neutrons, while the 
bottom two are for protons. In each panel the dotted, solid and dashed lines denote the distributions 
of the nucleons in the lsi/2, f'^id lpi/2 shell orbits, respectively. 
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FIG. 2. Same as Fig. [T| but for ^Li (left panels) and ^Li (right panels). 
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FIG. 3. Same as Fig. || but for ^^Li. The left and right panels are the results for the cases of 
S'n = 301 keV and 350 keV, respectively. 
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FIG. 4. Matter (pm), neutron (p„) and proton (pp) density distributions in ^"^■'^Li. The neutron- 
and proton density distributions arc obtained by summing up of the square of the corresponding 
single-particle wave functions. The matter distribution is simply given by pp + pn- 
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FIG. 5. /(y) for ^Li (left panels) and ''Li (right panels). In each column the top and bottom 
panels show those of neutrons and protons, respectively. In each panel the dotted, solid and dashed 
lines denote the distributions of the nucleons in lsi/2) lP3/2 and lpi/2 shell orbits, respectively. 
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FIG. 6. Same as Fig. |5| but for ^Li (left panels) and ^Li (right panels). 
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FIG. 7. Same as Fig. |5| but for ^^Li. The results shown in the left and right panels are for the 
cases of 5„ = 301 keV and 350 keV, respectively. 
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FIG. 8. Ratios of the structure functions of Li isotopes to that of the deuteron at = 10 
GeV^ without isoscalarity corrections. ^^Li-a(b) corresponds to the case of Sn = 301(350) keV. 



34 



2.0 



1.5 - 



1.0 



0.5 - 



0.0 



-1 I I ^ ^ I I ^ ^ r 



R/11-a) 

R^(l1-b) 




\ \ ^ 

■. \ \ 

\ \ ^ 
• \ 



\ \ ^ 
\ \ \ 

\ \ \ 



0.0 0.2 0.4 0.6 0.8 

X 



1.0 



FIG. 9. Rd{A) for Li isotopes (Q^ = lo GeV^). Rd{ll - a{b)) is for the case of ^^Li with 
Sn = 301(350) keV. 
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FIG. 10. Rn{A) for ^-^Li isotopes (Q^ = lo GeV^). The ratio of the difference between 
and F2 to F2 is also presented (denoted by Rn{D)). 
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FIG. 11. Rn{A) for ^^Li (Q^ = lo GeV^). i?„(ll - a(6)) is for the case of ^^Li with 
Sn = 301(350) keV. The ratio Rn{D) is also presented for comparison. 
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FIG. 14. Ratio of the multipole spin structure function, ( ^i, to the usual proton spin structure 
function, g^. The geometrical factor of 2\/5/3 is included in the denominator of the ratio. 
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